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Abstract 

The  present  study  generalizes  the  methodology  developed  in  [10]  to 
determine  visibility  probabilities  of  points  in  the  plane,  when  shadows 
are  cast  by  a  Poisson  field  of  random  objects  (disks),  given  a  source 
of  light  in  the  origin.  Furthermore,  a  technique  is  developed  for  the 
determination  of  the  moments  of  a  visibility  measure  on  star-shaped 
curves  in  the  plane.  A  numerical  example  is  provided  to  illustrate  the 
method  in  some  special  case. 

Key  Words:  i&JLb*JLL£y  Mediae,  Gzomz&U.caZ  P>wbabiLitie^ ,  Potion 
Random  Fi. nld.  Shadowing  Pkoczaa 


1.  Introduction 


In  a  previous  paper  [10]  we  studied  the  distribution  of  a  measure 
of  visibility  on  a  circle,  or  on  specified  arcs  on  the  circumference  of 
the  circle,  when  random  number  of  disks  of  random  size  are  randomly 
distributed  within  the  circle,  according  to  a  specified  Poisson  field. 
These  disks  obscure  the  visibility  (cast  shadows)  on  portions  of  the 
specif ieo  arcs,  for  an  observer  (source  of  light)  located  at  the  center. 

We  provided  a  method  of  computing  the  moments  of  the  measure  of  visibil¬ 
ity,  to  any  desired  order  by  a  recursive  technique.  We  also  approximated 
the  distribution  of  the  normalized  measure  of  visibility  on  (0,1)  by  a 
beta  distribution  mixed  with  a  two-point  discrete  distribution.  The 
present  study  provides  a  generalization  of  the  previous  methodology  to 
any  star-shaped  curves  in  the  plane,  when  the  random  disks  are  distri¬ 
buted  in  a  specified  region  according  to  a  Poisson  field.  Thi3  general¬ 
ization  can  treat  a  wide  class  of  geometrical  probability  problems  of 
line  of  sight,  visibility  measures  on  curves,  etc.  More  specifically, 
in  Section  2  we  introduce  the  general  structure  of  the  field  of  obstacles 
(disks)  and  the  associated  visibility  problems.  In  Section  3  we  further 
develop  the  basic  results  to  obtain  visibility  probabilities  of  points  on 
star-shaped  curves  in  the  plane.  Section  4  is  devoted  to  the  definition 
and  derivation  of  special  functions  (K-functions)  required  for  the  deter¬ 
mination  of  visibility  probabilities.  In  Section  5  we  discuss  the  measure 
of  visibility  on  smooth  star-shaped  curves.  We  develop,  furthermore,  re¬ 
cursive  formulae  for  the  computation  of  the  moments  of  the  measure  of 
visibility.  Section  6  provides  a  complete  theoretical  development  and 
numerical  computations  of  a  special  case  of  annular  region,  chosen  to 
illustrate  the  methodology  and  show  an  example  of  actual  computations. 

It  is  shown  that  the  mixed-beta  distribution  provides  an  excellent 
approximation  to  the  exact  distribution  of  the  measure  of  visibility,  in 
the  special  case  examined  in  Section  6. 

The  literature  on  shadowing  processes  is  quite  limited.  Chernoff 
and  Daly  [2]  studies  the  distribution  of  length  of  shadows  of  disks  on  a 
line.  The  shadowing  problem,  however,  is  a  special  case  of  the  general 
coverage  problem  on  which  there  is  extensive  literature.  In  particular 
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refer  to  the  studies  of  Robbins  [5,6],  Ailam  [1],  Greenberg  [4],  Siegel 
[7,8]  end  the  monograph  of  Solomon  [9]  which  summarizes  many  of  the 
important  results. 


Poisson  Random  Fields  and  Visibility  Probabilities 


Consider  a  countable  sec  of  disks,  V  ,  scattered  on  the  plane. 

Each  disk  is  specified  by  a  vector  (p,9,y)  ,  where  p  ,  0  S  p  <  »  , 
and  9  ,  -ir  £  9  <  ir  are  the  polar  coordinates  of  its  center  with 
respect  to  an  origin  0  (observation  point,  source  of  light);  and 
y  ,  0  s  y  <  y*  ,  is  its  diameter.  Let  denote  the  collection  of 

Bor el  subsets  of  the  space. 

SQ  *  {(P»9,y):  0£p<®,  -it£9 <ff,  0£y<y*}  . 

For  every  B  e  Sq  ,  let  N{B)  denote  the  number  of  disks  such  that 
(p>9,y)eB  (satisfying  condition  B).  If  N{B}  is  a  random  variable 
for  every  B  £  -Sq  then  the  elements  of  V  are  called  random  disks. 

In  particular,  the  elements  of  V  are  called  Poisson  random  disks  if 
for  every  B  e  Sq  ,  N{B}  has  a  Poisson  distribution  with  mean 

v{B}  -  y  /  /B/  dG  (y !  p ,  9 )  H  (dp ,  d9 )  , 

where  y  ,  0<y<®  ,  is  an  intensity  parameter;  H(p,9)  and  G(y|p,9) 
are,  respectively,  the  sigma-finite  measure  of  the  location  coordinates 
(p,9)  and  the  conditional  c.d.f.  of  the  diameter,  y  ,  given  (p,9)  . 

The  special  case  in  which  H(do,d9)  -  pdpde  and  G(y|p,9)  ■  G(y)  , 
is  called  the  standard  case.  The  standard  case  is  the  one  in  which  the 
centers  of  the  disks  are  uniformly  distributed  on  the  plane  and  their 
diameters  are  independent  of  their  location.  In  this  case 

y* 

v{B>  -  \if  Hg(y)dG(y)  , 

where  y  is  the  mean  number  of  disks  per  unit  area;  Hg(y)  is  the  area 
of  the  region  in  which  disks  of  diameter  y  ,  satisfying'  condition  B  , 
are  centered. 

The  theoretical  framework  is  not  restricted  to  random  disks  in  the 
plane,  but  can  be  generalized  to  Poisson  random  objects  in  the  plane, 
or  even  more  generally  in  the  space,  by  replacing  the  parametric  vectors 
(p,9,y)  with  a  general  parametric  vectors  (n^,...,np)  ,  which  charac¬ 
terize  the  location,  orientation,  shape  and  size  of  a  family  of  objects. 


In  such  cases  Sq  and  Sg  are  properly  generalized.  In  order  to  sim¬ 
plify  the  presentation,  we  restrict  attention  in  the  present  study  to 
the  family  of  random  disks  in  the  plane.  The  approach,  however,  is  valid 
for  general  families. 

A  natural  requirement  for  shadowing  processes  is  that  the  source  of 
light  (the  origin)  is  uncovered.  We  therefore  introduce  the  structural 
condition 

CQ  »  {(p,9,y);  ^£p<«,  -ir£0<ir,  0<y<y*}  (2.4) 

and  assume  that  all  random  disks  satisfy  Cq  .  In  special  cases  one  con¬ 
siders  more  stringent  structural  conditions,  specified  by  sets  C  con¬ 
tained  in  CQ  .  For  example,  in  the  previous  paper  of  Yadin  and  Zacks 
[10]  the  structural  condition  considered  is 

c  ■  {(p,9,y);  £p<l  -  -ffS0<ir,  0<y£l}  . 

Let  P  be  a  point  in  the  plane.  P  is  said  to  be  visible  (in  light) 
if  the  line  segment  OP  does  not  Intersect  any  random  disk.  The  set  of 
all  visible  points  in  direction  s  ,  -irss<ir  ,  starting  at  the  origin,  is 
called  a  line  of  sight,  L  .  Let  P  *  (r,s)  be  a  point  in  the  plane, 
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in  orientation  sand  distance  r  from  the  origin.  A  random  disk  (p,9,y) 
intersects  the  line  segment  OP  if,  and  only  if,  (p,9,y)  belongs  to 

B(r,s)  -  ((p,9,y);  (p,9)  e  B(r,s,y),  0<y£y*}  ,  (2.5) 

where  the  set  B(r,s,y)  is  the  set  of  all  points  having  distances  from 
0,P  smaller  than  y/2  .  (See  Fig.  1.)  A  disk  which  intersects  a  line 
segment  OP  is  said  to  cast  shadow  on  P  .  Accordingly,  a  point 
P  ■  (r,s)  is  visible  if  N(B(r,s)nC]  *  0  .  Hence,  a  line  of  sight  L 

s 

has  magnitude 

ML  II  ■  sup{r;  N{B(r,s)nC}  -  0}  .  (2.6) 

S 

Notice  that  ML  M  is  a  random  variable. 

3 

Let  P  ■  (r,s)  be  a  point  in  the  plane.  Under  the  Poisson  random¬ 
ness  assumption,  the  probability  that  P  is  visible  is 

Q(r,s)  -  P(N{B(r,s)nC}  -  0] 

-  exp(-v(B(r,s)nC}}  , 


(2.7) 


where  v{B(r,s)nC}  is  obtained  according  to  (2.2).  From  this  function 
one  can  immediately  obtain  the  distribution  of  I ! !  i  .  Indeed, 

P{ II L  1 1 >1)  -  Q(£,s),  0s^<»  . 

s 

Notice  that  in  the  standard  case 

v{B(r,s)nC0}  «  y£r  , 

where  £  is  the  expected  value  of  a  random  diameter,  Y  ,  whose  distri¬ 
bution  is  G(y)  .  Indeed,  for  Y  »  y  ,  the  area  of  Bq  *  B(r,s)nCQ  is 

H_  (y)  ■  ry  ,  as  can  be  implied  from  Fig.  1.  Accordingly,  the  distribu¬ 
te) 

tion  of  | |L  | |  in  the  standard  case  is  exponential  with  mean  1/p  , 
s 

where  p  ■_ u£  .  The  present  result  can  be  further  generalized  when  the 
random  objects  are  not  necessarily  disks  in  the  plane.  The  same  type 
of  result  for  infinitesimal  objects  is  given  in  Feller  [3,  p.10]. 

Let  P^  and  be  two  points  in  the  plane  having  coordinates 

(r(s^),s^)  and  (r(s2),S2)  •  A  necessary  and  sufficient  condition  for 
P^  and  P2  to  be  simultaneously  visible  is  that 

N{ (B(r(s1) ,s1)uB(r(s2))nC}  *  0  . 

Since  B(r(s^),s^)  and  B(r(s2),S2)  are  generally  not  disjoint  it  is 
more  convenient  to  consider  a  union  of  disjoint  sets  obtained  in  the 
following  manner.  Assume  that  the  points  (s,0,y)  in  C  have  ori¬ 
entation  coordinates  in  (s',s' Define,  for  any  s  in  [s',s'']  , 

B+(s)  -  {(p,9,y):  2sin(9-s')  <p<”»  3<9<s”j  0<y<y*} 

and 

B_(s)  -  {(p,9,y):  2sin^s-9 )'  <p<”»  s'<9<s,  0<y<y*}  . 

B+(s)  and  B_(s)  are  sets  of  all  points  in  C  ,  on  the  right  and  on  the 
left  of  the  ray  with  orientation  s  ,  which  do  not  intersect  it. 

Accordingly,  if  3^2  ,  (BCrCs^  ,s1)uB(r(s2)  ,s2>)nC  - 
C  -  (((B_(s1)uB+(s2))nC)u((B+(s^)nB_(s2))nC))  ,  where  A-B  is  the 
complement  of  B  with  respect  to  A.  Thus,  the  probability  that  the 
points  P^  and  P2  are  simultaneously  visible  is 

P(s.,s-)  -  exp{-[v{C}  -  v{(B  (s.)uB,  (s„))nC}  -  v{(BL(s.)nB  (s„))f!C}]}  . 
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Generally,  if  s '  £s^<. . .  <sq£s  '  *  are  the  orientation  coordinates  of  n 
points  then  the  probability  that  they  are  simultaneously  visible  is 

P(s^>»**,sq)  *  exp^^>{C}  -  v{  (B__ (s^)uB+(sq) )nC} 
n-1 

-  I  v{(B+(Si)nB_(si+1))nC}J|  . 


(2.13) 


Special  Field  Structure 

Consider  a  star-shaped  curve  in  the  plane,  C  ,  such  that  each  ray 
originating  at  the  origin,  0  ,  intersects  C  at  most  once.  The  curve 
C  is  specified  by  a  continuous  positive  function,  r(s)  ,  on  the  domain 
[s', s'*]  ,  where  -  y  £  s’  <  s”  s  ir/2  .  We  furthermore  require  that 
r(s)  will  have  almost  always  continuous  derivative. 

The  field  structure  is  specified  by  the  following  assumptions  on 
the  random  nature  of  the  distributions  of  the  disks.  First,  we  assume 
that  the  random  diameters,  Y^,Y2»...  of  the  disks  are  distributed  over  the 
interval  [a,b],  where  Osacbc*  .  Furthermore,  we  specify  two  star-shaped 
curves,  U  and  to  ,  between  the  origin,  0  ,  and  C  ,  so  that  the  cen¬ 
ters  of  the  disks  are  distributed  between  U  and  to  ,  and  no  disk  can 
either  cover  the  origin  or  intersect  C  .  More  specifically,  let 
u  -  (u(e);  e'sese”}  and  to  -  {w(9);  , 

where  [s'.s'1]  <=  [9,,0,r]  c  [-  y,  y]  .  We  impose  the  following  field 
structural  conditions  (see  Fig.  2):  Every  point  (p,0,y)  in  ,  with 
8  e  £9',9r ']  and  y  e  [a,b]  should  satisfy  the  two  conditions 

(C.l)  y  S  u(8)  <  p  £  w(8) 

(C.2)  If  a  point  P  *  (£(s),s)  is  on  the  circumf erence 
of  the  disk  l(£(s),s)  -  (w(9 ) , 9 ) i  ■  b/2  ,  then 

w(0)  cos  (9-s)  +  [(y)2  -  w2(0)  sin^CO-s)]1^2  <  r(s) 


Condition  (C.l)  guarantees  that  the  origin  cannot  be  covered.  Condition 
(C.2)  assures  that  C  is  not  intersected  by  any  random  disk  in  the 
field#  Since  -  y  £  s'  <  s''  £  ir/2  ,  B_(s')nC^  and  B+(s',)nC1  are 
disjoint.  Let  C  be  the  subset  of  C^  containing  only  disks  which 
may  cast  shadows  on  C  ,  i.e.,  they  intersect  line  segments  OP  ,  where 
PeC  .  According  to  (2.7),  (2.10)  and  (2.11)  ,  the  visibility  proba¬ 
bility  of  a  point  P  ■  (r(s),s)  on  is 

P^N{B(r(s),s)nC}  -  oj-  (3. 

exp£  (v{C}  -  v{B  U)nC}  -  v{BJ_(s)nC}]\  , 


Simultaneous  visibility  probabilities  of  several  points  on  C  are  given 
generally  by  formula  (2.13).  We  develop  now  the  technique  of  determining 
the  expected  number  of  points  in  the  sets  B  (s)nC  or  B  (s)nC 


4.  Definition  And  Derivation  of  the  K-Functions 

Given  a  point  P  »fc(s),s)  on  C  ,  let  yK_(s,t)  and  uK^(s,t) 
denote  the  expected  numbers  of  disks  ( p,8,y )  e  CL  with  s-tS0£s  and 
s58Ss+t  ,  respectively,  which  do  not  intersect  the  line  segment  OP  . 

These  functions  are  given  by 

s  w(8) 

K.(s,t)  -  /  /  G(y(p,s-0)|p,0)H(dp,d0)  , 

s-t  u(e ) 

and  (4.1) 

s+t  w(8) 

K  (s.t)  «  /  /  G(y(p,0-s) |p,8)H(dp,d0)  . 

s  u(8) 

where 

2psinr  ,  if  t<ir/2 

y(p.t)  »  (4.2) 

2p  ,  if  tiir/2 

Indeed,  when  1 9— s I  <ir/2  ,  a  disk  centered  at  (p,8)  does  not  inter¬ 
sect  OP  if  its  diameter  is  smaller  than  2psin|0-s|  .  If,  however, 
|8-s|Sir/2  ,  the  diameter  of  the  disk  should  not  exceed  2p  . 

Thus,  according  to  (3.2)  and  (4.1)  , 

v{C}  -  v{C1>  -  u[K_(s,,s'-8’)  +  K+(s,,,8”-s'')]  .  (4.3) 

Similarly,  for  all  s  e  [s',s'f]  , 

v{B  (s)nC>  -  y[K  (s,s-0’)  -  K  (s'.s'-e1)] 

(4.4) 

v{B+(s)nC}  -  y[K+(s,0"-s)  -  K+(s 1 '  ,0 '  '-s’ ' )  ]  . 

Furthermore,  for  s'  '  f  , 

v{B+(s1)nB__(s2)nC}  - 


S  —s  3  — s 

u£K+(sr  ”Vi)  +  k-(s2’  “V^)]  •  (4-5) 

S2~S1 

The  value  t  ■  — ^ —  was  selected  in  (4.5)  to  avoid  the  possibility  that 
a  disk  centered  at  (P,8)  with  s<0<s+t  ,  which  does  not  intersect  OP^ 
will  nevertheless  intersect  0P2  ,  and  vice  versa. 

For  the  purpose  of  evaluating  (4.1)  we  introduce  the  auxiliary 
functions 


and 
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v 

K  (s,d9,v)  *  /  G(y  (p , s-0 ) I p , 8)H(dp ,d0) 
0 


v 

K  (s,d0,v)  -  /  G(y(p,9-s)!p,9)H(dp,d8) 
0 


(4.6) 


The  K-functions  are  thus  given  by 

s 

K_(s,t)  -  /  (K_(s,d0,w(9))  -  K_(s,d0,u(9))] 

s-t 

and  (4.7) 

s+t 

K+(s,t)  -  /  [K+(s,d9,w(9))  -  K+(s,d8,u(9))]  . 

s 

Notice  that  G(y|p,9)  -  0  for  all  y<a  and  G(y|p,9)  -  1  for  all  y£b. 
Accordingly,  from  (4.6) 

A(0-s,v,b) 

K+(s,d6,\;)  -  /  G(y(p,0-s)  |  P,9)H(dP,d9)  + 

A(0-s,v,a) 


v 

f  H(dP,d0)  ,  (4.8) 

A(0-s,v,a) 

where 

°in(v’2ife7)  *  if  T<ir/2 

A(t,v,x)  •  (4.9) 

min(v,j)  ,  if  xS-rr/2 

The  function  K_(s,d9,v)  can  be  obtained  from  (4.8)  by  replacing  9-s 
by  s-0  .  In  the  standard  case  we  substitute  in  the  above  formula 
H(dp,d0)  ■  pdpd9  and  G(y|p,0)  ■  G(y)  .  As  an  example,  we  will  continue 
the  derivations  of  the  standard  case  and  a  uniform  distribution  of  Y 
over  [a,b]  ,  i.e., 

0  ,  y<a 

G(v)  -  .  a<v<h 
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Under  these  assumptions  one  obtains,  for  9>s  , 


A(9-s,v,b) 

K+(s,d0,v)  ■ j  (y(p,9-s)-a)pdp  + 

A(9-s,v,a) 


+  4p(v2  -  A2(8-s,v,b))  * 


0 


,  if 


9-s<sin 


:  3sin(9-s)  -  fv2 

* 

d0,  2  a2+ab+b2  , 

i  2  2  J 

!  12sin  (9-s) 

:  ^  r  2  a2+ab+b 2  ■, 

2  LV  12  J 


24  2 

z  sin (9-s) 


,  if 

,  if 

,  if 


sin_1(^7)S9-s<sin“1(^7) 
sin  "I'(^j)s9-s<Tr/2 
9-s*ir/2 


(4.11) 


As  seen  in  (4.11),  the  function  K^(s,d9,v)  ,  in  the  standard— uniform  case, 
be  written  as  K+(s,d0,v)  -  A(0-s, v,a,b)d9  .  Similarly,  K_(s,d9,v)  - 
A(s-9,v,a,b)d9  .where  A(9-s,v,a,b)  designates  the  R.d.S.  of  (4.11). 
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5.  The  Measure  of  Visibility  on  C  and  Its  Moments 

Define  the  indicator  function  I(s),  s'Ss^s''  so  that  I(s)  *  1 
if  (r(s),s)  is  a  visible  point  and  I(s)  -  0  otherwise.  The  total 
measure  of  the  visible  part  of  C  is  given  by 


V{C}  -  /  I(s)(r2(s)  +  (r’(s))2]  ds  . 

s' 


(5.1) 


V{C}  is  the  visibility  measure  of  C  .  In  coverage  problems  it  is  known 
also  as  the  measure  of  vacancy  (Ailam[l]>. 

The  first  moment  (expected  value)  of  V(C)  is 


E  (V(C) }  -  /  E{I(s)U(s)ds  , 


(5.2) 


2  r  2  * 

where  Z(s)  *  [r  (s)  +  (r  (a))  ]  .  Moreover, 


E{I(s) }  »  Q(r(s),s) 


expt-vtC^)  +  p[K_(s,s-0  )  + 


K+(s,9”-s)]} 


(5.3) 


V{C^}  is  given  by  (3.3)  In  the  standard  case  it  assumes  the  form 


[w2(0)  -  u2(9)]d9 

1  S' 


(5.4) 


Generally,  for  every  ns2  ,  the  n-th  moment  of  V{C}  is  given  by 


E{Va{C}}  -  /  ....  /  E{  II  I(s,)}  U  £(s  )ds 

s’  s'  1*1  1  i-1  11 


■  n  If  ...  f  P(s,,...,s  )  II  £(s.)ds. 

1  n  i-1  1  1 

t 

SSS.S.  .  .Ss  Ss’  ’ 

1  n 


(5.5) 


Indeed,  E{  II  I(s,)}  is  the  probability  that  all  the  n  points  are 
i-1  1 

simultaneously  visible.  According  to  (2.13),  (4.3),  (4.4),  and  (4.5)  we 


(5.6) 


6.  Numerical  Determination  of  the  Moments  of  Visibility  in  the  Standard- 
Uniform  Case  and  Annular  Sector  of  Disks 

In  Che  present  section  we  continue  the  development  of  the  formulae 
for  the  recursive  determination  of  the  moments,  when  the  disks  are  dis¬ 
tributed  over  the  annular  sector  with 


r(s)  -  r 
w(0)  ■  w 
u(9)  ■  u 


IT  |  ■  •  7T 

,  -  j  <  s'  s  s  s  s'  ’  <  J 


IT  TT 

,  -  yS  0'  S  6  S  0"  S  j 


ff  .  ,  .  IT 

,  -  y  S  0'  50  59”  S  y 


Here,  0<b/2<u<w<r-  b/2  . 

We  consider  the  standard  case,  where  the  distribution  of  the  disk  diam¬ 
eters,  Y  ,  is  uniform  (4.10).  Thus,  according  to  (4.7)  and  (4.11),  the 
K( s,t)  functions,  in  the  present  case,  do  not  depend  on  s  ,  and  satisfy 

/C+(s,t)  -  KJ s,t)  -  K*(t,w)  -  K*(t,u)  ,  ( 

for  all  s  ,  s ' Ssis ’ '  ,  where 


(6.1) 


K*(t,v)  *  C  A(r,v,a,b)dx 


in  which 


0 

,  if 

tSsin’1^) 

Kx(t,v) 

,  if 

sin’1  (~)  st<sin-1 

<2(t,v) 

,  if 

3in  1('^)st<?r/2 

^3(t,v) 

,  if 

ZZtt/2 

fv2  r 

'  2  2% 

1/2 

„  n  a  2  r 

•V  [ 

(4v  -a  ) 

• 

2v  cos  tl  -  y  v  It 

1-1 

■  ♦£[ 

(4v2-a 

2  1/2  1) 

,  )  -  a  cotan(t)  1 C 

(6.2) 


(6.3) 


K2CC,v)  »  -g 


-a  ^3  U 


(4v“-a“)  -  (4v“-b”) 


-  cotan  t)  , 

and 

<3(t,v)  -  <2(|,v)  +  (t  -  |)  [~-  -  -2^+fa2]  . 

Notice  that  in  the  present  case,  £(s)  »  r  for  all  s'ss^s''  .  Hence, 
the  moments  of  V{C}  can  be  determined  according  to  (5.4),  (5.7)-(5.9) 
and  (6.1)  in  the  following  manner. 


J-r  U 


sin  2TJ 


1  2  2 


X  -  exp{-  (w  -u)  (9* ’-9 ') >  , 


^(s)  -  exp{u[K*(s-0',w)  -  K*(s-0 '  ,u)  ]  }  , 


H(s)  •  exp{y[K*(0 ' *-s,w)  -  K* ( 0 f ' -s , u) ] } 


for  s'^s^s'' 

Define  recursively,  for  every  j2l  , 


9 

*j(a)  -  jt  4'j_1(y)exp{2u[K*(-2^  ,w)  -  K*(^  ,u)]}dy  . 


The  n-th  moment  of  V{C}  is  then 


u  *  \n!r  \  V  (s)H(s)ds  . 
n  4 ,  n-1 


Let  be  the  probability  that  V{C}  ■  L(C)  >  which  is  the  proba¬ 

bility  that  c  is  completely  visible.  According  to  (3.2)  and  (3.3) 

P2  ■  X  exp{jj(K*(s’-0 '  ,w)  -  K*(s'-0',u)  + 

K*(0"-sM,w)  -  K*(0"-s",u)]> 
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The  numerical  determination  of  the  moments  u  was  performed  in  the 

EX 

following  maimer.  The  interval  [s', s'']  was  partitioned  to  M  sub¬ 
intervals  of  length  A  ■  (s''-s')/H  .  The  grid  points  for  the  computa¬ 
tion  are  then  Sj«s'+jA  ,  j»0,l,...,M  .  The  functions  ¥q(s)  end  G(s) 
were  computed  exactly  for  each  .  The  functions  4'k(s),kil  , 

were  then  computed  at  the  grid  points  s.  ,  according  to  the  numerical 


integration  formula 


i 


to  .if  j-0 

2  (Vl(so)exp{2v  *  (2)}  +  Vi(sl))  ’  if  j-1  (6>12) 

W  *  ')  2  <Vl(so)exp{2w  *  (^)}  +  Vl(3l))  + 

+  A  jE  ¥.  1(s.)exp{2u  K  (i±llA) }  ,  if  ji2 

i-1  1  1 

where  ^(t)  ■  K*(t,w)  -  K*(t,u)  .  Finally,  the  n-th  moment  of  V(C)  is 
determined  by  the  formula 


Xn!rnA{(¥  .(s  )H(s  )  +  ¥  .  (s  )H(s  ))/2 

n-i  o  o  n— l  m  m 


mrl 

+  2  ¥  i (s,)H  (s .) }  .  (6.13) 

j-1  a“1  j  i 

If  the  value  of  M  is  large  one  obtains  very  good  approximation  by 

applying  formulae  (6.12)-(6.13) .  In  Table  6.1  we  provide  the  first 

ten  normalized  moments  of  V(C)  ,  i.e.,  u  /(r(s"-s'))n  ,  for  the  case 

n 

of  s'  ■  -ir/18,  s''  ■  -s'  ,  M*60  ,  r»l  ,  u*.5  ,  w».75  ,  a«.l  and  b*.3  . 

The  values  of  u  (the  Poisson  intensity)  are  1(2)9  . 
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Table  6.1  The  Normalized  Moments  of  V(C)  and 


Approximations ,  for  s'»"7l8,  s" 


w-».75.  a«.l,  b*.3 


Their  Beta-Mixture 


M«60,  r-1. 


Intensity 


1 

2 

3 

4 

Normalized  Moments 

5  6  7 

8 

9 

10 

00 

0.951 

0.934 

0.926 

0.921 

0.918 

0.917 

0.917 

0.919 

0.921 

0.925 

.901 

0.951 

0.934 

0.926 

0.921 

0.917 

0.915 

0.913 

0.912 

0.910 

0.910 

0.860 

0.816 

0.794 

0.781 

01773 

0.768 

0.765 ' 

0;763- 

0.763 

0.765 

.730 

0.860 

0.816 

0.794 

0.780 

0.772 

0.765 

0.761 

0.757 

0.754 

0.752 

0.77C 

0.713 

0.681 

0.662 

0.650 

0.643 

0.638 

0.634 

0.633 

0.633 

.592 

0.778 

0.713 

0.681 

0.662 

0.649 

0.640 

0.634 

0.629 

0.625 

0.621 

0.703 

0.623 

0.584 

0.562 

0.548 

0.538 

0.532 

0.527 

0.525 

0.524 

.480 

0.703 

0.623 

0.584 

0.562 

0.547 

0.536 

0.528 

0.522 

0.518 

0.514 

0.636 

0.545 

0.502 

0.477 

0.461 

0.451 

0.444 

0.439 

0.435 

0.434 

.389 

0.636 

0.545 

0.502 

0.477 

0.460 

0.449 

0.440 

0.434 

0.429 

0.425 

Notice  that  u  «P,  . 

00  1 

As  in  the  previous  study  of  Yadia  and  Zacks  [10],  the  distribution  of 
V{C}/r(sf '-s')  is  approximated  by  a  mixture  of  a  beta-distribution  with 
a  discrete  distribution  concentrated  on  0  and  1  .  The  mixed  distri¬ 
bution  of  the  normalized  visibility  measure,  V  {C}  »  V{C}/r(s' ’-s’)  , 
has  a  c.d.f. 

r0  ,  if  x<0 

,  x 


F*  (x) 


V(l-Vpi>  5 (S 


■*  o 

hr l 


ya_1(l-y)3”1dy  ,  if  0sx<l 


,  if  x^l 


(6.14) 


Since  the  value  of  is  known,  we  determine  f  (x)  ,  in  the  various 

cases,  by  equating  the  first  three  moments  of  (6.14)  to  the  moments  of 
* 

V  {Q  •  The  n-th  moment  of  (6.14)  is 


P.  +  (1-Pn-P.)  7 


q(q+l) . . .  (a-fn-l) 


0  V  (q+S)  (q+0+1) . . . (q+B+n-1) 


,  n*l 


(6.15) 


Accordingly,  if  cQ  *  un“?]_  ,  ttien 


2C;- 


clc3-c2 


<  Vc2>  (c2-c3) 


C1C3“C2 


wr 


c1(ct+6) 


if  if  y  ' 

In  Table  6.2  we  provide  the  values  of  a  =■  ) 


,  a  and  3 


corresponding  to  the  cases  of  Table  6.1.  We  have  also  computed  the 
first  ten  moments  of  (6.14)  according  to  the  parameters  determined  by 

A 

the  true  moments,  u  ,  and  presented  them  in  Table  6.1  (lower  line  of 

Q 

each  case) . 

We  see  that  the  actual  moments  and  the  ones  of  the  approximating  dis¬ 
tribution  are  extremely  close.  This  indicates  that  the  mixed-beta 
distribution  is  apparently  very  close  to  the  true  one.  We  obtain  in 
this  manner  also  an  estimate  of  the  probability  of  complete  coverage, 
PQ  .  An  explicit  formula  for  this  parameter  is  not  yet  available. 

Table  6.2  Parameters  of  the  Mixed-Beta  Approximating  Distribution 
(Specifying  parameters  as  in  Table  6.1) 

Intensity  ""  - 

u  o  Pn  P,  o  8 


0 

po 

P1 

a 

8 

.1733 

.0020 

.9005 

1.1161 

1.0398 

.2765 

.0152 

.7302 

1.1405 

1.0967 

.3290 

.0326 

.5921 

1.0956 

1.1194 

.3590 

.0536 

.4801 

1.0450 

1.1386 

.3754 

.0768 

.3893 

0.9937 

1.0967 

* 


‘4 


r 

t 


K 


t 


In  order  Co  realize  Che  efface  of  differenc  case  paramecers  on  Che 

resules,  we  presenc  in  Table  6.3  che  first  10  momencs,  u  ,  of 
*  ® 

V  (C)  ,  for  a  longer  curve  (s'»-it/3,  s'^-s')  .  The  compucacion 

was  performed  on  a  grid  of  H^ISO  subincervals.  The  incensicies 

W  are  1(2)9,  and  che  discribucion  of  che  disk  diameCers 

is  uniform  on  (0.1, 0.3). 

We  do  noc  presenc  Che  momencs  of  (6.14),  which  correspond  Co  chese 
cases,  since  chey  are  very  close  Co  che  accual  momencs.  The  resules 
of  che  compucacions  presenCed  in  Table  6.3  show  similar  crends  co 
chose  indicaced  in  Tables  6.1  and  6. .2. 

* 

Table  6.3  normalized  Momencs,  V  (C)  ,  And  The  Paramecers  of  The  Mixed 
Beca  Discribucion  for  che  Cases  Of  s'»~ir/3.  s,,*-s',  r»l. , 
ix«.5,  w».75,  a».l»  b*.3  and  M»180. 


_1L_ 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

OB 

1 

.951 

.912 

.880 

.854 

.832 

.814 

.799 

.787 

.776 

.767 

.686 

3 

.860 

.758 

.682 

.623 

;577 

.541 

.511 

.487 

.467 

.451 

.322 

5 

.778 

.631 

.529 

.456 

.402 

.360 

.328 

.303 

.282 

.266 

.152 

7 

.703 

.526 

.412 

.335 

.280 

.241 

.212 

.189 

.171 

.157 

.071 

9 

.636 

.438 

.321 

.246 

.196 

.162 

.137 

.113 

.104 

.094 

.033 

Table  6.5  Paramecers  of  The  Mixed-Beca  ApproximaCing  Discribucion 
(Specifying  Paramecers  as  in  Table  6.3) 


u 

0 

po 

P1 

3 

1 

.0856 

0 

.6856 

15.8793 

2.9360 

3 

.1364 

0 

.3223 

8.9669 

2.3286 

5 

.1619 

0 

.1515 

6.3166 

2.2383 

7 

.1762 

0 

.0712 

4.9572 

2.3287 

9 

.1838 

0 

.0335 

4.1394 

2.5105 

r 

(■ 
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